A compact star rotating at 1122 Hz and the r-mode instability 
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ABSTRACT 

We show that r-mode instabilities severely constraint the composition of a 
compact star rotating at a sub-millisecond period. In particular, the only vi- 
able astrophysical scenario for such an object, present inside the Low Mass X- 
ray Binary associated with the x-ray transient XTE J 1739-285, is that it has 
a strangeness content. Since previous analysis indicate that hyperonic stars or 
stars containing a kaon condensate are not good candidates, the only remaining 
possibility is that such an object is either a strange quark star or a hybrid quark- 
hadron star. We also discuss under which conditions sub-millisecond pulsars are 
rare. 

Subject headings: star: neutron — star: rotation — r-mode — quark star - 
submillisecond pulsar — gravitational waves 



Very recently iKaaret et al.l (120071 1 reported the evidence of a X-ray transient with a 
pulsed component of the emission having a frequency / = 11 22 ±0.3 Hz. This signal is 
interpreted as due to the rotation of the central neutron star. As such this object would be 
the most rapidly rotating compact star discovered up to now. This single observation clearly 
needs to be confirmed, maybe by the analysis of future X-ray transients of the same object. 
The implications of this rapid rotation on the Equation of State (E OS) and in particular on 
the allowed values of the m ass and the radius have been discussed in lLavagetto et al.l (120061 ) 



and in lBejger et al.l (120071 ). Here we discuss the problem of t he stability respect to r-modes of 



such a rapidly rotating object (for a review on r-modes see lAndersson fc Kokkotad (120011 )) 



It has been widely discussed in the literature the possibility that r-mode instabilities 
can very efficiently drag angular momentum from a rotating compact star if its temporal 
evolution in the fl — T plane ( angul ar ve locity and temperature) ente rs the r-mode instability 
window, see e.g. lAndersson I (119981 ) and lFriedman fc Morsinkl (119981 ). Therefore huge regions 
of the Q — T plane are excluded. Moreover, the size and position of that window is strictly 
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related to the composition of the star, since it is strongly dependent on the value of bulk and 
shear viscosity. It is particularly important to recall that for stars containing strangeness, 
as hyperon stars (Lindblom fc Owen! 120021 ) . hybrid stars (IDrago et al.l 120051 ) and strange 
quark stars (IMadsenl 120001 ) there is also a contribution to the bulk viscosity associated with 
the formation of strangeness. Due to this, the instability window splits into two parts: 
one which starts at temperatures larger than (7 ± 3) 10 9 K (High Temperature Instability 
Window HTIW) and a lower temperature window at temperatures smaller than (5 ±4) 10 8 K 
(Low Temperature Instability Window, LTIW) 0- Concerning the left border of the LTIW, 
its position is regulated by the shear viscosity and by the so called viscous boundary layers 
locat ed at the interface between t he crust and the fluid composing the inner part of the 
star (IBildsten &: Ushomirskyl 120001 ) . In particular, in bare quark stars, either composed by 
normal or by superconducting quark matter, due to the absence of a significant crust, the 
left border of the LTIW extends to much lower temperatures and the LTIW has a minimum 
corresponding to a significantly lower temperature than in the case of stars containing a 
crust. 

In this Letter we discuss in which region of the Q — T plane the compact star originating 
XTE J 1739-285 is most likely located, due to its composition. Let us start by discussing 
the simplest possibility, i.e. that the object is a neutron star. In principle a neutron star 
can rapidly rotate in two cases, either if its temperature is very large, above a few MeV, 
or if it is recycling, spinning up due to mass accretion (see Fig. 1). Concerning the first 
possibility, a hot neutron star would be a newly born one, since the time needed to cool 
below one MeV is of the order of one minute. This is clearly not the case of the stellar 
object under discussion. Concerning recycling, it should take place on the left side of the 
in stability window, l ocate d at lower temperatures (see Fi g. 1). The main result of the analysis 
of I Andersson et al.l (120001 ). revisiting previous analysis (lLevinlll999l ; IBildsten fc Ushomirsky 
20001 ) is that a neutron star can never spin- up to a rotational period shorter than ~ 1.5 ms. 
This result is based on the estimates of the temperature and of the mass accretion rates 
of Low Mass X-ray Binaries (IBrown fc Bildsten Ill998l ; iBhattacharyyal 120021 ) indicating that 
temperatures lower than ~ 10 8 K cannot be reached. Therefore a sub-millisecond neutron 
star cannot be present at the center of a Low Mass X-ray Binary (LMXB). It is important 
to remark that this conclusion is confirmed by more recent analysis, taking i nto account the 



composition stratification of the rigid crust (IGlampedakis fc Andersson 



2006) and discussing 



the nonlinear development of the r-mode instability (IBondarescu et al.l 120071 ). 

Let us now discuss the case in which the compact star contains strangeness. It has been 
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shown, in th e case of hyperon stars (jWagonerl 120021 ; iReisenegger fc Bonacidl2003l ). and for 
hybrid stars (IDrago et al.ll2006l ) (strange quark stars also have a similar behavior), that due 
to the large value of the bulk viscosity associated with the non-leptonic weak decays, there 
are two windows of instability, the LTI W and the HTIW introduced above. The HTIW does 
not affect significantly the angular velocity of the star because the cooling of a newly born 
star is so fast that there is not enough time for the r-mode instability to drag a significant 
fraction of the angular momentum. Therefore, the star exits the HTIW with an angular 
velocity close to the initial one. When the temperature drops down to a few 10 s K the star 
reaches the LTIW and it starts to lose angular momentum due to r-mode instability. In 
Fig. 1 we show examples of the instability windows in the case of a pure quark star and of a 
hybrid star. Notice that the position of the LTIW depends rather strongly on the mass of the 
strange quark m s . For large values of m s the instability windows shrink considerably. In our 
analysis we have considered two possibilities concerning the value of m s : a small value ~ 100 
MeV, of the order of the strange quark current mass, and a large value ~ 300 MeV, similar 
to what has been obtained in NJL-like models of quar k matter in which a density-dependent 
constituent mass has been computed (IBuballa II2005I ). 



Concerning the left side of the LTIW, it is easy to see that it cannot be used to accom- 
modate a submillisecond pulsar. Indeed, in the case of hybrid or hyperonic stars the left side 
of the LTIW is similar to the one of neutron stars, discussed above, and we can apply here 
the same analysis done for a neutron star and concerning the minimal temperature of the 
compact object at the center of a LMXB. Therefore also hybrid or hyperonic stars cannot 
recycle to frequencies significantly exceeding 700-800 Hz if they are to the left of the LTIW. 
Moreover, in the case of quark stars the absence of viscous boundary layers implies that they 
can rotate rapidly only on the right side of the LTIW (a very light crust suspended over 
the surface of a quark star is irrelevant from the viewpoint of r-mode instabili t y, alth ough 
it can be important for the production of x-ray bursts, see iPage fc Cummingl (120051 )). In 
conclusion, stars containing strangeness can rotate at submillisecond periods only if they are 
to the right of the LTIW. 

Let us now discuss the results of our analysis also con sidering the constraint s posed 
onto the EO S by th e mass shedding limit. The main result of lLavagetto et al.l (120061 ) and of 
Bejger et al.l (120071 ) is that soft EOSs, and in particular the ones based on hyperonic matter 
or on matter in which kaon condensation takes place, are rather unlikely. In fact, not only 
the range of allowed masses is rather small, but moreover the only confi gurations satisfyin g 
the stability constraint for mass shedding turn out to be supra massive (IBejger et al.l 120071 ). 
making it very difficult to use these stellar structures in a mass accretion scenario. This resul t 
confirms previous discussions of other astrophysical objects ruling out soft EOSs (IOzelll2006l ). 
On the other hand, stars containing deconfined quark matter are not excluded if the quark 
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EOS is stiff enough (lAlford et al.l 120071 ). We can therefore conclude that the object under 
discussion is either a quark or a hybrid star. 

We can now study the temporal evolutions of the angular velocity of the star Q, of 
its temperature T and of the amplitude of the r-modes a whic h are calculated by solving 
the set of differential equations given in lAndersson et al.l (120021 ) (Eqs. 15-23-24). The only 
technical difference in our calculation is the inclusion of the reheating asso ciated wi t h the 



dissipation of r-modes by bulk viscosity, as d iscussed in various papers (jWagonerl 12002 



Reisenegger Bonacid 120031 : iDrago et al.ll2006l ). Due to the reheating, the time needed for 
the star to cool down along the right border of the LT IW is rather long. In the equation 
regulating the thermal evolution (lAndersson et al.ll2002l ) we have added a term: 



E, 



bv 



2E/r bv 



X 



where E ~ 10 51 a 2 Mi.4i? 2 Pr 3 2 erg is the energy of the r-modes (the mass, radius and ro- 
tation period of the star are here expressed in units of 1.4 M Q , 10 km and in millisecond, 
respectively) and r bv is the bulk viscosity. Due to the large value of the bulk viscosity for 
strange matter at temperatures of a few 10 8 K, the related dissipation of r-modes strongly 
reheats the star. The trajectory in the Q — T plane describing the time evolution of the star 
follows essentially the border of the LTIW and the star keeps rotating as a sub-millisecond 
pulsar for a very long time, strongly dependent on the value of m s and ranging from ~ 10 2 
years when m s = 100 MeV to ~ 10 8 years when m s = 300 MeV, as shown in Fig. 2. In Fig. 3 
we show the temporal dependence of the r-mode amplitude. Notice that, in the absence of 
reheating the time needed to slow-down to frequencies below 1 kHz is at maximum of the 
order of ~ 10 4 years, and can be much shorter than a year if a small value of m s is adopted. 
In the past this branch of the instability window was only discussed in connection with very 
young pulsars because reheating due to bulk viscosity was not taken into account. 

Before discussing the possible astrophysical scenarios for a submillisecond pulsar inside 
a LMXB, it is important to remark two constraints that a realistic model should fulfil. First, 
a LMBX is an old object, with a typical age of 10 8 -10 9 years which is the time needed by the 
companion (having a low mass) to fill its Roche lobe and to start accreting on the neutron 
star. Therefore we need to provide a mechanism allowing the central object inside the LMXB 
to rotate rapidly while being so old. Second, submillisecond pulsars are certainly rather rare, 
and in particular there is no evidence of a uniform distribution of pulsar rotational frequencies 
extending from a few hundred Hz up to more than a kHz. Therefore a realistic model should 
also indicate why most of the compact stars will not be detected rotating at a submillisecond 
period. As we will show, these two constraints can be satisfied in two different astrophysical 
scenarios: a first scenario in which a old hybrid or quark star is accelerated up to frequencies 
exceeding one kHz by mass accretion and a second scenario in which the quark or hybrid 
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star is born with a submillisecond period and it is now spinning down by r-mode instability 
(see Fig. 4). 

The possible realization of these scenarios depends on two main ingredients: 

• the value of m s , which, as we have already shown, regulates the magnitude and position 
of the LTIW; 

• the cooling rate of the central region of the star which determines the inner temperature 
T of a star which is accreting material from the companion (IMiralda-Escude et al. 



199(1 ) . 



Concerning the cooling rate, if strange quark matter is present in a compact star, direct 
URCA processes are possibl e and therefore the cooling is (generally) fast. It turns out from 
Miralda-Escude et al.l (119901 ) that in this case the inner temperature is ~ 5 x 10 7 K for a mass 
accretion rate of M — lCr lo M0/year and that the temperature scales as T oc M 1//6 . Another 
possibility recently proposed is that, due to the formation of diquar k condensates, URCA 
processes are strongly suppressed and the cooling turns out to be slow (IBlaschke et al.ll2000l ). 
At the same time, bulk viscos ity can still be large enough to suppress r-mode instability 
(IBlaschke &: Berdermannl 120071 ) . In the following we discuss both possibilities, either of a 
fast or of a slow cooling. 

Let us consider first the case in which m s = 300 Me V. 

The first scenario, i.e. the star span up by accretion, can indeed be realized if T > 10 8 K 
(see lower panel of Fig. 1). Such a temperature can be reached via reheating due to mass 
accretion. In the case of fast cooling a rather large value of mass accretion rate is needed, 
M ~ lO~ 8 M /year, and this stringent request can explain why submillisecond pulsars are 
rare. A model in which the cooling is slow is instead excluded, because it would be extremely 
easy to re-accelerate the star to very large frequencies and therefore submillisecond pulsars 
would not be rare. 

Also the second scenario in which the star is spinning down due to r-mode instabilities is 
possible. As shown in Fig. 2, the time spent by the star above 1122 Hz is of order of 10 8 years 
and it is therefore compatible with the typical age of LMXBs. In this case submillisecond 
pulsars are rare because only a (small) fraction of ne wly born compact stars can r otate 
wi th submillise c ond p eriods (see e.g. the discussions in lAndersson fc Kokkotasl (120011 ) and 



in 



Perna et al.l (120071 )). In this scenario the distribution of millisecond and submillisecond 



pulsars is determined by the mass of the star at birth: if the mass is large enough, strange 
quark matter can form in the core of the star, r-modes are efficiently damped by the large 
bulk viscosity and the star can be a submillisecond pulsar at birth; if the mass of the star 
is small, strange quark matter can not form, the bulk viscosity of nucleonic matter can not 
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damp r-modes and the refore in a short period the star loses almost all its angular momentum 



flLindblom et al]ll998h . 



We consider now the case in which m s = 100 MeV. 

The first scenario in which the star is spinning up by accretion can be realized but T must be 
at least ~ 5 x 10 8 K for the star to be accelerated up to 1122 (see upper panel of Fig. 4). Such 
a temperature can be reached but it needs strong reheating due to a large mass accretion 
rate and slow cooling. 

The second scenario, i.e. the star is spinning down by r-mode instability, is ruled out since, 
as shown in Fig. 2, the time spent by the star above 1122 Hz is of order of 100 years and 
this is incompatible with the typical age of a LMXB. 

In conclusion, the two astrophysical scenarios proposed to explain a submillisecond 
pulsar inside a LMXB can both be realized by a quark or a hybrid star if m s ~ 300 MeV. 
Instead, if the m s is small, only the first scenario is possible and it requires a slow cooling. 
Our results are summarized in Table 1. 

The main uncertainties in our analysis are due to the possible existence of other damping 
mechanisms, taking place on the left side of the LTIW. For ins tance magnetic fields can be 



important to suppress r-mode instabilities (IRezzolla et al.ll2000i ). but their effect is probably 



negligible for frequencies exceeding ~ 0.35 Qk, if the internal magnetic field is not larger than 
~ 10 16 G. Obviously, even larger uncertainties exist concerning quark matter. As discussed 
above, the bulk viscosity of quark matter strongly d epends on the strange q uark mass and 



on the possible formation of a diquark condensate (lAlford fc Schmittl 120071 ). Clearly, our 



analysis can provide much needed constraints on the EOS of quark matter. 

Finally, let us stress that the outcome of our analysis is that a compact star rotating at 
a submillisecond period inside a LMXB can only be a quark or a hybrid star. Future obser- 
vations will be important to clarify if the object at the center of XTE J1739-285 constitutes 
indeed an example. 



It is a pleasure to thanks Carmine Cuofano for many useful discussions. 
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Table 1: Microphysics and astrophysical scenarios 



m s [MeV] 


cooling 


scenarios (I or II) 


conclusions 


300 


fast 


I: O.K., M > 10- 8 M Q /y 

II: O.K., needs rapid rotation at birth 


O.K. for 
both scenarios 


300 


slow 


I: no, submillisecond pulsars too frequent 
II: O.K. 


excluded 

due to scenario I 


100 


fast 


I: no, T is too low 

II: no, slow down too rapid 


not possible 

in both scenarios 


100 


slow- 


I: O.K., needs large M 

II: no, slow down too rapid 


O.K. for 
scenario I 



Table 2: Connection between microphysics (mass of the strange quark and cooling rate) 
and astrophysical scenarios for a submillisecond pulsar inside a LMXB. The first scenario 
corresponds to a compact star which is spinning up due to mass accretion; the second scenario 
to a compact star born as a submillisecond pulsar and still in the process of slowing down. 
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Fig. 1. — r-mode instability windows. The star is unstable and loses angular momentum 
by emitting gravitational waves in the regions above the instability lines displayed in figure. 
Thin solid lines correspond to neutron stars , for t wo extreme values of masses and radii 
allowed by the analysis done by iBej ger et al.l (120071 ). The shaded area in between is repre- 
sentative of intermediate values of masses and radii. Dashed lines delimit the LTIW and 
the HTIW for strange quark stars. Thick solid lines correspond to hybrid stars. In both 
cases the mass and radius of the star are M = 2M and R = 13 km. In the upper panel 
m s = 100 MeV and in the lower panel m.„ = 300 MeV. The bulk and shear viscosities are 
in general from lAndersson fe Kokkotad f 20011 ) . Bulk viscosity of strange quark matter at 



high temperatures is from iMadsenl J200ol )~ The arrow indicates the minimal temperature 
considered possible in LMXBs. The dot indicates the actual position of the compact star 
associated with XTE J1739-285, as it results from our analysis. 
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Fig. 2. — Time dependence of the rotational frequency of the compact star. The thin lines 
take into account only the reheating due to shear viscosity, the thick lines take into account 
also the effect of reheating due to bulk viscosity. We show results for two different values 
of the strange quark mass, m s =100 MeV and m s =300 MeV. The dashed horizontal line 
corresponds to a frequency of 1122 Hz. 
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a 




Fig. 3. — Temporal evolution of the amplitude a of the r-modes. Results 
without reheating in panel (a) and taking into account reheating in panel 
(b). The amplitude of the r-mode is significantly damped when reheating is 
taken into account, the maximum amplitude is reduced by roughly 2 orders of 
magnitude. Therefore non linear effects should be suppressed. 
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Fig. 4. — Upper panel: first scenario (see text). The star is born with a low angular velocity. 
After a cooling period at a constant angular velocity (line (a)), the star reaches the instability 
window and loses angular momentum during a long period, ~ 10 9 years (line (b)). After 
this long period the mass accretion process starts and the star is heated (line (c)) and re- 
accelerated (line (d)) to a frequency of 1122 Hz (indicated by the thick dot). Lower panel: 
second scenario (see text). The star is born with a large angular velocity. After a short 
period of cooling at a constant frequency (line (a)) the star enters the instability window 
and loses angular momentum via r-mode (line (b)). The frequency of the star can be larger 
than 1122 Hz (indicated by the thick dot) for a period of ~ 10 8 years. The mass and radius 
of the star are M = 1.4M and R = 10 km, respectively. 



